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Abstract. We investigate the thermal emission from magnetic neutron star surfaces in which the cohesive effects of the magnetic 
field have produced the condensation of the atmosphere and the external layers. This may happen for sufficiently cool (T < 
10 6 ) atmospheres with moderately intense magnetic fields (about 10 13 G for Fe atmospheres). The thermal emission from an 
isothermal bare surface of a neutron star shows no remarkable spectral features, but it is significantly depressed at energies 
below some threshold energy. However, since the thermal conductivity is very different in the normal and parallel directions to 
the magnetic field lines, the presence of the magnetic field is expected to produce a highly anisotropic temperature distribution, 
depending on the magnetic field geometry. In this case, the observed flux of such an object looks very similar to a BB spectrum, 
but depressed by a nearly constant factor at all energies. This results in a systematic underestimation of the area of the emitter 
(and therefore its size) by a factor 5-10 (2-3). 
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1. Introduction 

Almost a decade after the discovery of the soft X-ray source 
RXJ 1 85635-3754 fWalter, Wolk, & Neuhauser 1996) using the 
ROSAT-PSPC, the thermal component associated with the 
direct emission from a neutron star's surface has been de- 
tected in more than 20 X-ray sources. In many cases it is su- 
perimposed on a power-law tail, but seven of these objects 
are well-characterized as simple blackbodies with tempera- 
tures ranging between 60 and 100 eV. The apparent small 
emitting surface of RXJ 185635-3754, inferred from the best 
blackbody fit and its parallax ( Walter & Lattimer 2002 1 have 
led to speculation about its nature, and whether its appar- 
ent smallness can be considered as evidence that the object 
is a strange star: a self-bound object made of up, down and 
strange quarks ( Po ns et al. 2 002 Dra ke et al. 2002> . or a stan- 
dard, misaligned pulsar. Although this new observational class 
(isolated compact stars) is probably the first real opportu- 
nity to place stringent constraints on the equation of state 
(EOS) of dense matter from astrophysical measurements (see 
ILattimer' & Prakash 2001 for a review), one must be cautious 
before concluding that an apparently small X-ray source is a 
quark star, because the X-ray spectrum is not the only infor- 
mation available. Using the Hubble Space Telescope, Walter 
& Matthews ( 1997 > subsequently identified an optical source 
at 6060 A and 3000 A, with a brightness only about 7 times 
brighter than an extrapolation of a 62 eV X-ray blackbody 
into the optical V band. The observed optical fluxes have 
been confirmed by subsequent observations from the 2-meter 
NTT dNeuhatiser, Thomas, & Walter 1998| and new HST mea- 
surements (Pons et al. 2002 1. Remarkably, the other three iso- 



lated compact X-ray sources that have been detected in the 
optical band (RX J0720.4-3125, RX J1308.6+2127, and RX 
J 1605. 3 +3249) also have a significant optical excess over the 
extrapolation of the X-ray blackbody (a factor 5 to 14). None 
of them have yet been detected as 1 .4 GHz radio sources. Thus, 
RXJ 185635-3754 is not an uncommon object, for it shares the 
same general observational properties of other isolated neu- 
tron stars (blackbody spectrum in X-ray, no evident spectral 
features, optical excess), except for the fact that four of them 
have fairly long periods (8-22 s), while RXJ 185635-3754 is not 
variable, with a reported upper limit on the pulse amplitude of 
< 1.3% (IBurwitz et al. 20031 

Since blackbodies are no more than a simple approxima- 
tion of the true emission mechanism, a step forward in under- 
standing the thermal emission of neutron stars, consists of com- 
puting model atmospheres for low magnetic fields (< 10 n G) 
and emergent spectra for several compositions, masses and ra- 
dius (iRorriani 19871 IMillerT992l |Rajagopal & Romani 1996) . 
The parameters that determine the shape of the observed spec- 
trum are the atmospheric composition, the red-shifted temper- 
ature, To, = r e ff/(l + z), the redshift z, the interstellar medium 
absorption, iih, and the ratio Rco/d, where d is the distance 
to the object and R^ - R(l + z), R being the radius of the 
star. For a Planck spectrum the redshift factors contribute only 
to an overall scale factor, so that the redshift does not affect 
the models. For realistic atmospheres with spectral features, 
however, the redshift is an important parameter because the 
identification of spectral lines would determine the M/R ratio 
and establish relevant constraints on the EOS of dense mat- 
ter. After a detailed investigation of the different parameters 
that affect the observed spectra, Pons et al. (2002) concluded 
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that X-ray data alone do not allow us to establish severe con- 
straints, given the large number of degrees of freedom. The 
combination of the X-ray spectra with the optical fluxes turned 
out to be much more restrictive, as was suggested earlier by 
|Pavlov et al. (1996*)| The broadband spectral energy distribu- 
tion (SED) of isothermal heavy element atmospheres can fit the 
observed SED, but these models predict detectable absorption 
lines and edges. It was hoped that the more sensitive spectro- 
graphs on Chandra and XMM-Newton would be able to de- 
tect spectral features and help to understand the nature of this 
objects, but no lines or absorption edges have been seen in 
long exposures (IBurwitz et al. 20011 IDrake et al. 20 02 >. Only 
the presence of broadband departures from a single blackbody 
spectrum has recently been suggested ( Bur witz et al. 2 003 1. 

The puzzle, then, is as follows: single temperature models 
based on heavy element atmospheres explain the broadband 
spectrum and the lack of variability of the X-ray spectrum, 
but they are in contradiction with the absence of spectral fea- 
tures. Single blackbody models, besides being unrealistic, can- 
not account for the systematic optical excess that is observed 
in 4 sources. An anisotropic surface temperature distribution 
(the simplest case would be two blackbody components) can 
reconcile X-ray and optical observations, because the quantity 
(Roo/D) 2 must be interpreted as the solid angle subtended by 
the star's surface area which is visible at some distance D. Thus 
the true value of will generally be larger if the assump- 
tion that the temperature is uniform on the surface is relaxed. 
However, this is barely consistent with the lack of pulsation 
of RXJ 185635-3754, unless its magnetic and rotation axes are 
oriented in a very particular way. 

Nevertheless, there is an important ingredient missing in 
the above discussion, the likely existence of strong magnetic 
fields. The strong (10 I2 ~ 13 G) magnetic fields of pulsars suggest 
that most, if not all, neutron stars should have similarly strong 
fields. However, detailed atmosphere models with strong mag- 
netic fields are only available for hydrogen (Zav lin et al. 19 95 1, 
because reliable opacities and EOS have not yet been de- 
veloped for heavier elements and because of the complexity 
of modeling magnetic atmospheres with arbitrary magnetic 
field structures. For heavy-element dominated atmospheres, 
only approximate treatments of magnetic Fe atmospheres ex- 
ist (Rajagopal, Romani, & Miller 1997 1, and the results show 
that the spectra are globally much closer to a blackbody than 
for light element atmospheres. There is a case to be made that 
RX J 185635-3754 is magnetized. The deep VLT image of the 
target released by van Kerkwijk & Kulkarni ( 2000 1 shows what 
looks like a classic bow-shock nebula. The presence of a bow- 
shock suggests that this is a magnetized neutron star with a rel- 
ativistic wind, as seen in the pulsars PSR 1957+20 (Kulkarni & 
Hester lT988l and PSR 2224+65 (Cordes, Romani, & Lundgren 
fT993l . 

An alternative possibility that explains naturally the ab- 
sence of spectral features is the emission from a solid surface. 
This was a common idea 20-30 years ago JBrinkmannT980l 
B80 in the following) until the existence of a thin gaseous 
atmosphere was appreciated and model atmospheres became 
more popular. However, at sufficiently low temperatures, 
highly magnetized neutron stars may undergo a phase transi- 



tion that turns the gaseous atmosphere into a solid (La i 200 ft . 
The critical temperature below which the atmosphere conden- 
sates depends on the composition and the magnetic field. For 
example, for typical magnetic field strengths of 10 13 G, a Fe 
atmosphere will condensate for T < 0.1 keV while a H atmo- 
sphere needs temperatures lower than 0.03 keV to undergo the 
phase transition to the metallic state (Lai 2001 1. Notice that ef- 
fective temperatures of the observed isolated neutron stars fall 
in this temperature range, therefore they should plausibly be 
in the solid state if the dominant element in the atmosphere is 
Fe. In such a metallic neutron star surface made of nuclei with 
atomic number Z and atomic weight A, the pressure vanishes 
at a finite density 

p s * 560 AZ- 3/5 B^ 5 g cirT 3 (1) 

where B\2 is the magnetic field in units of 10 12 G. 

This idea has been recently revisited by Turolla et al. 
(2004), who found that the emitted spectrum is strongly de- 
pressed at low energies, thus making more difficult the recon- 
ciliation between observational data and the condensed surface 
model. However, in this last work the emissivity is calculated 
neglecting one of the transmission modes in some cases. This 
simplification is inaccurate and can modify the emitted spec- 
trum, as we will discuss in the text. Besides, the sole presence 
of a strong magnetic field, necessary to condensate the atmo- 
sphere, results naturally in an anisotropic surface temperature 
( Geppert, Kiiker, & Page 2004 ), which must be calculated con- 
sistently with the magnetic field structure. In this paper our aim 
is to study the thermal emission from solid surfaces of neutron 
stars and its implications on the observational properties of iso- 
lated compact objects. 

The paper is organized as follows: Sect. 2 is devoted to dis- 
cussing the calculation of the emissivity starting from the de- 
scription of the dielectric tensor for the condensed neutron star 
surface. In Sect. 3 we show results, including integrated spec- 
tra for different strengths and geometries of the magnetic field. 
Final remarks and a summary of the main conclusions are dis- 
cussed in Sect. 4. 

2. Equations and input microphysics 

2.1. Dielectric Tensor 

The emission properties (reflectivity, emissivity) of a neutron 
star surface in a condensed state are obtained from the dielec- 
tric tensor, and were first studied in detail in B80. In the pres- 
ence of the strong magnetic fields expected to exist at the sur- 
face of a neutron star, the dielectric tensor changes significantly 
with respect to the non-magnetic case, leading to birefringence 
and other associated phenomena. In order to calculate the di- 
electric tensor in the presence of magnetic fields and to take 
into account dissipative processes, it is better to obtain first the 
conductivity tensor cry, and then to calculate the dielectric ten- 
sor €ij according to the equation 

An 

eij = 5ij + i — a- u . (2) 

CO 

In this expression, to is the angular frequency of the electro- 
magnetic wave propagating in the medium. 
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For a magnetized plasma, electrical conductivities for the 
static case (u> = 0) have been calculated by some authors 
( |Hernquist 1984| IPotekhin"T^ 99 1 with and without taking into 
account the quantizing effect of the magnetic field. The general 
expression for the conductivity tensor can be written as follows 



CTi 



f 2 Mb ( df \ 



(3) 



where e is the electron energy, / the Fermi-Dirac distribution, 
Nb is given by the expression (see Potekhin 1999 for details) 

"max 

N B = ' . f-, /. g n V(e/c) 2 - (m e c) 2 - 2m e Haj B n, (4) 

and the functions ry play the role of effective relaxation times. 
For a non-quantizing magnetic field, these effective relaxation 
times can be expressed in terms of a relaxation time, to, which 
is the inverse of the collisional frequency v D (D after damping) 
However, in the quantizing case, collisional frequencies in the 
longitudinal and perpendicular direction to the magnetic 
field are no longer equal. For a quantizing magnetic field, the 
components of the Ty tensor are 



i + (cj b t ± ) 2 



Tyx 



1 + {CJ B T±) 2 



(5) 



with T|| and Tj_ being the inverse of the effective collisional fre- 
quencies, vj^ and vj, respectively, and o>b = ^7 being the elec- 
tron cyclotron frequency. 

In the dynamic case (a> + 0), the conductivity tensor can be 
obtained by the transformation 



T, - 10) 



T11 - ItO 



(6) 
(7) 



and the dielectric tensor, for degenerate non-relativistic elec- 
trons, reads 



4tt 

■■ 6 i} + i—a-ij 



( S -iD ' 
iD S 
OP, 



(8) 



where we have defined the following complex quantities 

,2 



R 



= 1 - 



(jJ- 



11) 



P = 1 - 



OJ U) + Ll>b 

2 

art 



ID 2 + !WV? 



S = -(R + L); 
D = \{R~L). 



(9) 

(10) 

(11) 
(12) 



Here, u> p = (4jre 2 n e /m e ) 1 ^ 2 is the electron plasma frequency 
and, using the density at zero pressure given by Eq. Q, the 
electron particle density n e can be calculated as follows 



n e = 1.24 x 10 27 Z^ 5 B% 5 cm" 3 . 



(13) 



2.2. The dispersion relation. 

Consider an electromagnetic wave propagating in a medium 
described by a dielectric tensor e a B- The dispersion relation is 
easily obtained by introducing the Maxwell tensor 



2 w 
Wj = kikj — k 5ij + 



(14) 



where k\ are the Cartesian components of the wave vector. 
Since, in terms of the Maxwell or dispersion tensor, the electric 
field of the wave satisfies the equation 



KjEj = 0, 



(15) 



the condition to have a non-trivial solution leads to the disper- 
sion equation 



Ay [= 0. 



(16) 



In the same Cartesian frame as used in B80, where the mag- 
netic field B is in the x — z plane and forms an angle a with 
the z-axis, which is normal to the surface, the dielectric tensor, 
obtained by a rotation transformation from Eq. 0, adopts the 
form 



' S cos 2 a + P sin 2 a -iD cos a sin a cos a(P - S ) 

iD cos or S -iD sin a 

k sin a cos a(P - S ) iD sin a P cos 2 a + S sin 2 a 



(17) 



For our purposes it is convenient to write the dispersion rela- 
tion for the transmitted wave in terms of the incident one. We 
consider an incident wave with wave vector 



ki = kj(- sin i cos sin i sin /3, - cos z), 



(18) 



where i is the angle of incidence and /3 the azimuth. Since the 
wave is propagating in vacuum and we are neglecting vacuum 
polarization, the dispersion relation gives k, = at/c. The re- 
flected and transmitted waves have wave-vectors 

k r = ki(- sin i cos sin i sin j3, cos i) 

k, = k,{- sin6 m cos/3, sin6>,„ sin/3, - cos# m ); m = 1,2, (19) 

where the subscript m refers to the ordinary and extraordinary 
modes and m is the angle of refraction, which is a complex 
number. Introducing the complex refractive index for the trans- 
mitted wave n = kc/co, and using the Snell law, sin^ = 
sin ;'/«„„ the dispersion relation leads to a quartic equation for 
the refractive index n, 

n 4 (P + v sin 2 a) + n 2 (gv - 2PS + u sin 2 a) + PRL + gu = 

sin i sin(2a) cosy3(n 2 - sin 2 i) l/2 (u + n 2 v) (20) 

where v — S — P, u — PS — RL, and g - sin 2 z'[l - sin 2 a(l + 
cos 2 /?)]. This quartic equation can be solved analytically, after 
squaring, and the refractive indexes for both ordinary and ex- 
traordinary modes can be obtained after analyzing which are 
the two physical roots that satisfy the original Eq. (1201 . 

In Fig.[^we show the real and imaginary parts of the two 
modes as a function of the photon energy and for two different 
magnetic fields B = 5 x 10 12 G (solid lines) and fi = 5x 10° 
G (dashes). The corresponding electron plasma frequencies are 
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Fig. 1. Refractive indexes without (left) and with (right) damping effects. The real parts are shown in the top panels and the 
imaginary parts are shown at the bottom. Two different values of the magnetic field are depicted, 5 x 10 12 G(solid lines) and 
5 x 10 13 G (dashed lines). The photon incidence angle has been taken to 43° to compare with Fig. 6 of Turolla et al. 2004 



3.6 keV and 14.3 keV, respectively. We have taken a magnetic 
field normal to the emitting surface and an incidence angle of 
43°. The left panels correspond to the case without damping 
(yf = V _L = 0) while the right panels show the equivalent re- 
sults including effects of collisional damping (these collisional 
damping frequencies were calculated by using the public code 
developed by Potekhin '). The resonance at the plasma fre- 
quency u>p is clearly visible, as well as a second, narrow reso- 
nance at u>b (~ 58 keV for B - 5 x 10 12 G). More interestingly, 
we can observe a region at low energy in which one of the 
modes takes a large imaginary part, which will result in sig- 
nificant absorption. The energy at which this happens can be 
estimated as follows. For simplicity, let us assume normal in- 
cidence and a magnetic field normal to the surface, and let us 
consider the case without damping. In this limit the dispersion 
relation is 

Pn 4 - IPS n 2 + PRL = 0, 
and the roots are simply n 2 = L and n 2 — R , i.e., 



0){<jJ ± 0)b) 

1 www. iof f e . rssi . ru/astro/conduct/condmag . html 



Therefore, for energies u> < ^ [-1 + -y/1 + 4(w p /w B ) 2 J the 
second mode acquires an increasingly larger imaginary part. 
It turns out that for B > 10 12 G we are always in the case 
a>p <sc o>b, and the above condition becomes approximately 
a> < (tPjJuB- In brief, we can expect two main features in the 
spectrum, a resonant absorption near the plasma and cyclotron 
frequencies, and lower emission (compared to the BB) below a 
certain threshold energy * oj 2 /(l>b- In Turolla et al. (2004i, the 
modes with imaginary part larger than a certain value (0.01) 
were neglected, this is indicated in the figure by the horizontal 
dashed line. In some cases, this may result in lower emissivity 
at low frequencies, and we prefer not to neglect them. 

2.3. Reflectivity and emissivity 

Knowing the complex refractive indices of the two modes, and 
solving Eq. d!5i . the ratio between the electric field compo- 
nents, (E' x , EL EL), of the transmitted wave can be obtained: 

Y± = a m = [D sin a + S sin i — n m sin i cos B 
+ (P cos 2 a + S sin 2 a)(« 2 - S - sin 2 i sin 2 B)\ x 
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Fig. 2. Emissivity normalized to the BB emission as a function 
of energy for B = 10 12 G and T = 10 6 K and varying the 
incident angle. The magnetic field has been taken to be normal 
to the emitting surface. 
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Fig. 3. Normalized emissivity integrated over all possible inci- 
dent angles as a function of energy forfi = 10 12 GandT = 10 6 
K. Different orientations of the magnetic field are compared: 
from top to bottom = 4.5, 18, 36, 45, 54, 63, 72; and 85.5°. 



sin i cos) 



+ iD sin i sin B (cos a sjn'^- 



sin i + sin a cos a(P - S ) 



(S ~ n 2 m ) + 



sin a cos a 



sin i + sin a sin i cos B | + 

-l 



and 



= b„ 



[(P-S)sin 2 /sin 2 y8 + D 2 ]) 
|a m (sin 2 i cos B sin/? - iD cos a) 



(22) 



sin B sin i ^jn^ - sin 2 i + iD sin a 
x (sin 2 B sin 2 i - n 2 + S ) . 



(23) 



These relative amplitudes coincide, after some minor algebraic 
manipulation, with those in B80. In Turolla et al. (2004i, a dif- 
ferent linear combination of the three equations arising from 
Eq. (II 51 was used, but both results can be shown to be equiva- 
lent when the refraction index is a solution of Eq. (I20> . 

The boundary conditions at the surface of separation of 
both media, (vacuum and the magnetized plasma) imply that 
the tangential component of the magnetic and electric field and 
the normal component of the magnetic and the electric dis- 
placement must be continuous. This results in the following 
equation 



E ± 

E\\ 
E» 



Si 



B, 



ft] COS 6\B\ »2 cos ^2^2 

cos i cos i 



\ 


(E'n) 




E ' 2z 




El 


/ 


v E » ) 



(24) 



where E, E' and E" are, respectively, the electric field of the 
incident, transmitted and reflected wave, and the subscripts ||, 
_L refer to components parallel or perpendicular to the inci- 
dence plane. Above, we have defined A m = b m sin B - a m cos B, 



B m = b m cosB + a m smB, and C m = e 3 ia m + e$ib m + e 33 . This 
system of equations can be solved for the electric field of the 
reflected wave in terms of the incident one. Details and explicit 
expressions are given in the Appendix. 

The reflectivity can be now obtained by assuming that the 
incident wave is the incoherent sum of two linearly polarized 
waves parallel and perpendicular to the incidence plane. From 
Eq. (IA.U and taking E ± = 1, E\\ = 0, the two complex compo- 
nents of the reflected field (E^ and E ± ) can be calculated, and 
the reflectivity is simply 



p ± = |£J 2 + |<| 2 . 



(25) 



Analogously, for the case of polarization in the incidence plane, 
one takes E ± = 0, E\\ - 1, and after obtaining E,, and E ± the 
reflectivity is given by 



P\\ = \E\ 



| 2 + |<| 2 . 



(26) 



To obtain the reflectivity of the unpolarized incident, p v , wave 
we just take the average 



Pv = ^(Px+P\\) 



(27) 



Finally, the reflectivity is related to the emitted specific in- 
tensity I v by Kirchhoff's law 



/ v = (l-p v )B v , 

where B v is the Planck intensity 

2hv 3 /c 2 



B v = 



(28) 



(29) 



e\p(hv/kT) - 1 

In Fig. |2] we illustrate how the normalized emissivity, i.e., 
emitted intensity normalized to the BB value, a v = I Y /B V = 
1 - p v , varies with the angle of incidence for a fixed magnetic 
field of B = 10 12 G normal to the surface. The temperature is 
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Fig. 4. Same as Fig.|3]but for B = 10 13 G. 

T = 10 6 K. The emissivity is strongly reduced compared to the 
blackbody case for energies lower than 0.2-0.5 keV, depending 
on the incident angle. A reduction for energies close to that cor- 
responding to the plasma frequency (indicated by the vertical 
line) is also evident in all cases. Notice how for incident angles 
close to 7r/2 the emission is strongly suppressed in a wide range 
of energies. 

Let us now comment on the variation of the emissivity 
with other relevant parameters. In Fig. we show the normal- 
ized emissivity, integrated over all possible incident angles, as 
a function of the photon energy. We have taken T = 10 6 K 
and B = 10 12 G and we show results for different magnetic 
field orientations (a is the angle between the magnetic field 
and a vector normal to the surface). At energies greater than 
the electron plasma frequency (about 1 keV in this case) the 
emitted flux approaches the BB value (a v = 1), but the spec- 
trum is significantly depressed at low energies. As explained 
above, this is due to the fact that the refractive index has a large 
imaginary part. A resonance, which produces a reduction of 
the emissivity, is also visible at energies close to the electron 
plasma frequency. Notice that the resonant energy is not ex- 
actly the plasma frequency but depends on the angle between 
the magnetic field and the surface. This can be understood con- 
sidering that the refractive index of one of the modes becomes 
infinity when the coefficient of n 4 in the dispersion relation 
vanishes. Again, neglecting damping for simplicity, we have 
P + (S - P) sin 2 a = 0, which leads to 



co 4 - (co B + d) 2 p )(i) 2 + co 2 co 2 B cos 2 a = 0. 
The solution of this biquartic equation is 



u>\ + w 2 



1 + 



1 



4(l> 2 . o) 2 B cos 2 a ^ 
p p 

(oj 2 b + w 2 ) 2 



1/2 



(30) 



(31) 



For magnetic fields perpendicular to the surface a — 0, we have 




0.001 0.01 0.1 

E (KeV) 

Fig. 5. Same as Fig.|4]but including the effect of the motion of 
ions. The vertical line corresponds to the ion cyclotron energy. 

plasma frequency in the a = case is shifted to lower and 
lower energies. 

The same but for B = 10 13 G is shown in Fig. |4] 
Qualitatively we obtain the same behaviour, but the plasma fre- 
quency is larger and the resonant energy is shifted consistently 
to higher values. Notice that the frequency below which the 
spectrum is depressed depends weakly on the magnetic field 
u> 2 /cl>b oc B 1 ^ , and for standard values of the magnetic field 
falls in the range 0.1-0.2 keV. In this energy range the effect of 
interstellar medium absorption can make it difficult to distin- 
guish between the two effects. 



2.4. The effect of motion of ions. 

When we were finishing this paper, our attention was drawn to 
a very recent preprint (van Adelsberg, Lai, & Potekhin 2004 1 
with similar results to ours. In this work, the authors included 
terms related to the motion of ions in the dielectric tensor, 
which result in smaller reflectivity (larger emissivity) at fre- 
quencies below the ion cyclotron frequency = — ). The 
way in which the effect of ions is included is a crude simplifi- 
cation: as if they were free ions although they are in a lattice. 
It is quite doubtful that this approximation can actually rep- 
resent reality, but it gives some interesting results that are an 
indication that more work is needed along that line. In short, 
free ions can be introduced in the calculations just by modify- 
ing the components of the dielectric tensor (|8) in the following 
way ( |Ginzburg 1970||van Adelsberg, Lai, & Potekhin 2004| i: 



1 



(to + (jl>b){Ci> + (l>Bi) + iur/^A* 

7 1 

cot + OJZ, 



D " 



(32) 
(33) 



in 2 and oj 2 



but for magnetic fields parallel to the where a> p . = (47rZ 2 e 2 n,/m i ) I/2 is the ion plasma frequency and 



emitting surface one finds a> 2 = w 2 +oj 2 and oj 2 - 0. Therefore, 
as we approach a - n/2, the resonance corresponding to the 



U>B i nie 

A = 1 + — (1 -Z _1 )+ — . 



(34) 
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In Fig. |3 we show the normalized emissivity integrated over 
all incident angles, for T = 10 6 K and B = 10 13 G, including 
the effect of ions. The vertical line indicates the energy corre- 
sponding to the ion cyclotron frequency. 

In order to understand the effect of ions we follow the same 
argument after Eq. (12 1 i . Neglecting damping for simplicity, 
and considering that the ion plasma frequency is much smaller 
than the electron plasma frequency (a> Pj <k u> p ), the dispersion 
relation is simply 



\R ) (0J + O>b)(C0 + OJBi) 

When to > there is not much difference with respect to the 
case in which ions are not included. But when u < cl>b,, the new 
term in the denominator changes sign and the refractive index 
of the mode that acquired a large imaginary part when the ion 
contribution was not included becomes 

n 2 = 1 + p —, (36) 

and the mode is no longer damped. 



1000 




0.001 0.01 0.1 1 

E (keV) 

Fig. 6. Integrated emitted flux, L v /4nD 2 , (in arbitrary units) for 
a uniform temperature of T = 10 b K. The magnetic field ge- 
ometry has been taken to be a dipolar distribution and sev- 
eral values of the magnetic field strength are compared: B p = 
10 12 ,5xl0 12 , 10 13 , and 5xl0 13 G (values at the magnetic pole). 



3. Spectral energy distribution. 

The total spectral emission of the star, the specific luminosity, 
is obtained by integrating the specific intensity over the solid 
angle pointing outward and over the surface of the star, 

r*7.7i r*n fin 

L V = R 2 I d(f> I sin GdO I d<p' (37) 
Jo Jo Jo 

dff /„(£/, (p, 0', (j)') cos 0' sin 0' 

Jo 

where and <p are the polar and azimuthal angles of a given 
point at the surface and 0' and 4>' the polar and azimuthal angles 
of the direction defining the solid angle element. 

The presence of the magnetic field, however, makes the 
emission anisotropic, and the observed specific flux F Y is differ- 
ent from the integrated emitted flux obtained from the specific 
luminosity, i.e. L v /4nd 2 for a star at a distance d. To calculate 
the observed specific flux we have to integrate only over the 
observed hemisphere (OH). Denoting by R the radius of the 
star, 

R 2 C 

F v = — \ I v {6, cos & sin6» d0d<f> (38) 

d Joh 

where the angles ff and <p' specify the angular direction from 
the surface element towards the observer. Note that in this pa- 
per we are neglecting general relativistic effects. For example, 
light bending increases the area of the observed hemisphere, 
depending on the compactness (M/R) of the neutron star, as 
discussed in Page d!995> or Psaltis, Ozel, & DeDeo ©00 1. 

In order to get a more accurate feeling of what a realistic 
emission would be, in the following we will consider that the 
magnetic field has a dipolar geometry, this is 

B(R, 6) = y (1 + 3 cos 2 0) 1/2 (39) 



with Bp being the intensity of the magnetic field at the pole. 

The integrated spectrum as a function of the energy is 
shown in Fig.|6]for different intensities of the dipolar magnetic 
field and a temperature of T = 10 6 K. For reference, the BB 
spectrum is also depicted (solid line). The broadband spectrum 
is essentially featureless and the flux is systematically lower 
than the BB flux for the same temperature; this effect is more 
significant at lower energies. Notice that this flux must be un- 
derstood as an average over the whole emitting surface, which 
is radiating anisotropically. The real observed flux depends on 
the particular location of the observer, as defined by Eq. (1381 . 
In Fig. we show the flux observed from three different di- 
rections, forming an angle of 0°, 45°, and 90°, with the polar 
axis, respectively. Near the maximum the differences between 
observers can be at most a factor 2, but the observer location 
becomes less relevant as we move towards either the low and 
high energy part of the spectrum. 

At this point, the two main features of the "metallic sur- 
face" model seem to be: first, an almost featureless spec- 
trum, and second, an overall flux smaller than that of a 
BB at the same temperature, especially at low energies. 
However, as we mentioned in the introduction, there is an 
important point to consider. The fact that the condensed sur- 
face is strongly magnetized makes the thermal conductivity 
very different in the directions perpendicular and parallel to 
the magnetic field lines. Similar effects have been pointed 
out to be relevant in the envelope (Greenstein & Hartke 
1983 Page 1995 1 or in the crust, where a very recent study 
(Geppert, Kuker, & Page 2004 1, finds that the anisotropy in the 
temperature distribution depends very strongly on the partic- 
ular geometry of the internal magnetic field, resulting in vari- 
ations of temperature of up to a factor 5. In a separate work 
( Perez- Azorin, Miralles & Pons 2 005 1, we will report results 
from a detailed study of the temperature distribution obtained 
from 2D diffusion calculations for different magnetic field ge- 
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Fig. 7. Observed flux, F v , (in arbitrary units) for a uniform tem- 
perature of T = 10 6 K and three different observation angles. 
The magnetic field geometry has been taken to be a dipolar 
distribution with B p - 5x 10 13 G. 



ometries. For the purpose of understanding qualitatively the ef- 
fects on the observed spectrum, in this paper we will limit our 
analysis to the case in which the temperature distribution has 
the following angular dependence 



T = T p [cos 2 6 B +xsm 2 6 B ] 



1/4 



(40) 



where 9b is the angle between the field and the normal to the 
surface, x is me rat i° between the thermal conductivities nor- 
mal and parallel to the magnetic field, and T p is the polar tem- 
perature (where 9b = 0). The origin of this distribution has 
been discussed in previous works on neutron star envelopes 
(Greenstein & Hartke l 19831 Page ll995> . Note that for a dipolar 
magnetic field ^ is a function of the polar angle because the 
magnetic field strength varies with the latitude. 

In Fig. [8] we show the same three cases as in Fig.^but for 
the anisotropic temperature distribution given by Eq. J40> with 
T p = 10 6 K. The optical band of the spectrum is not very much 
altered, but the high energy tail is significantly depressed. This 
effect, combined with the low energy depression caused by the 
high reflectivity of the metallic surface at low energies, results 
in a broadband spectrum that mimics the BB spectrum but with 
an overall reduced flux by nearly a factor of 10. This means 
that, for a fixed distance to the source, the observed flux from 
such a particular neutron star surface will look like a Planckian 
spectrum, but the apparent area of the source (and therefore the 
radius) would be underestimated by a large factor. To make this 
point more explicit, in Fig|5]we plot the observed flux (dashes) 
from a model with B p = 5 x 10 13 , T p = 10 6 and 9„ = 90°, as 
seen after taking into account interstellar medium absorption 
with n H — 1.4 x 10 2() cm -2 , compared with a uniform temper- 
ature, blackbody model that fits the X-ray part of the spectrum 
(solid line). The parameters of the BB model are: T = 10 6 
K, Hh = 1.3 x 10 20 crrT 2 , and a relative normalization fac- 
tor (oc (Rcc/d) 2 ) of 1/5. Consequently, the apparent estimated 
value of the R^/d is 2.3 times lower than that of the "real" 
model, despite the X-ray spectrum being very similar. 
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Fig. 8. Observed flux, F v , (in arbitrary units) for an anisotropic 
temperature distribution described in Eq. J40l > with T p — 10 6 
K and three different observation angles. The magnetic field 
geometry has been taken to be a dipole with B p = 5 x 10 13 G. 
The corresponding Planckian spectrum for T = 10 6 K is also 
shown for comparison. 



For comparison, we have also included the effects of ions in 
Fig- 13 (dotted line). Notice that the optical flux of this model is 
a factor 3 larger than the BB fit prediction, similarly to what 
has been observed in isolated neutron stars. Nevertheless, it 
must be stressed that this is just a crude approximation, with 
no pretension of being the real answer to the observed optical 
excess, but it serves to illustrate how important it is to under- 
stand details about the magnetic field structure, the properties 
of the solid lattice, and the temperature distribution, before one 
is able to make robust estimates of the neutron star properties 
(e.g. radius). 

4. Final remarks. 

In this paper we have revisited the bare neutron star sur- 
face model first studied in detail in B80, which in the 
last years is becoming popular again and attracting the 
attention of other groups (Lai 2001 Turolla et al. 2004 
van Adelsberg, Lai, & Potekhin 2004). Our results for constant 
temperature magnetized surfaces confirm qualitatively those 
reported earlier by other authors, with small quantitative dif- 
ferences due to the fact that in previous works (Turolla et 
al. |2004 1 some approximations (neglecting one mode) were 
made, which made the results dependent on the cutoff value 
of the imaginary part of the refraction index. In general, mod- 
els with uniform temperature show a broadband spectrum that 
is very close to Planckian at energies above « 2 /wb, and signif- 
icantly depressed (up to a factor 10) in the optical band. The 
spectrum is almost featureless, with only some small bumps 
at energies where the interstellar medium absorption makes 
it difficult to distinguis and fine tune between different pa- 
rameters. However, in our opinion, there is a key point that 
is barely addressed in previous works and needs more atten- 
tion: in the crust and the condensed outer layer the assump- 
tion of a homogeneous temperature distribution is inconsistent 
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Fig. 9. Observed flux, F v , (in arbitrary units) for the anisotropic 
temperature distribution described in Eq. (I40i with T p = 10 6 K 
and O = 90°, as it would be seen after taking into account inter- 
stellar medium absorption with tin = 1.4 x 10 20 cirT 2 . We show 
results with (dots) and without (dashes) including the effect of 
ions. A uniform temperature, blackbody fit of the X-ray part 
of the spectrum is also depicted with solid lines (T = 10 6 K, 
rin- 1-3 x 10 20 crrT 2 ) but corrected by a factor 1 /5, that is, the 
apparent estimated value of the R^jd is 2.23 times lower than 
that of the "real" model. 

because magnetic fields of the order of or larger than 10 13 G 
imply some degree of anisotropy in the thermal conductivity. 
The analogous effect in the inner crust has recently been stud- 
ied ( |Geppert, Kiiker, & P age 2004), finding large variations of 
temperature when the magnetic field is confined to the crust. 
Notice that, if in the inner regions we have superconducting 
protons, as seems to be the case, this situation is very likely. A 
detailed analysis of the realistic, self-consistent emission from 
bare neutron star surfaces requires, therefore, multidimensional 
transport calculations with the presence of magnetic fields, and 
using appropriate boundary conditions (accordingly with the 
calculated a v ). Such calculations, as well as fits to real data, are 
in progress ( P erez-Azorin, M iralles & Pons 2005 ) and will be 
reported elsewhere. In the mean time, one can guess what sort 
of changes to expect by looking at the emitted spectrum pro- 
duced by an ad-hoc temperature distribution, as we discussed 
in this paper. This example was very illustrative of one fact: 
the observed flux of such an object is very close to a BB spec- 
trum, but we might be underestimating the area of the emitter 
(and therefore its size) by a large factor. In addition, depending 
on the strength of the magnetic field, and including the effects 
of ions, we could even obtain an optical flux larger (relative to 
the BB case) than that in the X-ray band, which is commonly 
found in all isolated neutron stars with an optical counterpart. 

As stated before, general relativistic effects have not been 
included in this work. A first correction is simply to redshift 
all energies and temperatures. This is of crucial importance if 
spectral features are present, but only translates into an overall 
scale factor if we consider a BB or a featureless spectrum. A 
second effect might be more relevant. Light bending increases 
the observed emitting area, smearing out partially the differ- 



ences between different observers, depending on how compact 
the object is. In order to make precise parameter estimates at 
least these two major corrections should be included. In sum- 
mary, there is much physics to be understood and analyzed in 
the near future before drawing robust conclusions on the nature 
of isolated compact objects (neutron stars vs. strange stars), and 
measuring with precision their radii and masses. 
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Appendix A: Reflected wave amplitudes. 

The system of equations d24l> can be solved for the electric field 
of the reflected wave in terms of the incident one. The result 
reads as follows: 

E" ± = D l {[Bi(l - wi)(A 2 sin i + C 2 cos i)- 

B 2 (l - Wz) (Ai sin i + C\ cos ;')] E ± + 
2B\B 2 {w\ -W2)sinz'cos2'E|[} 

E[ = £T> {2(A 1 C 2 -A 2 C 1 )E ± + 

+ [B\(l + w\) (A 2 sin i - C 2 cos i) - 

B 2 (l + w 2 )(Ai sin z-Ci cos z)]£||} (A.l) 

where 



V n ».- Sin 1 

w m = : , (A.2) 

COS I 

D = Bi(l + w\) (A 2 sin/ + C 2 cos i) - 

B 2 (l + w 2 ) (Ai sin i + C\ cos i) (A.3) 

From these equations the original expressions in B80 can be 
recovered by using the relation 

C m - sin 2 i + A m w m sin i cos i. (A. 4) 

We have found a typo in the first line of Eq. (19) in B80, where 
the last J{ 2 should be a S 2 . Notice that Turolla et al. (2004i 
expressed the amplitudes of the reflected waves in a different 
way. The reason for this apparent difference is the same that 
led to different expressions for the coefficient a m defined by 
Eq. (I22> : they used a different linear combination of the three 
equations arising from Eq. (I15> . We have checked numerically 
that both formulations are equivalent when n m is a root of Eq. 
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